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The phase structure of the two-flavor Polyakov-loop extended Nambu-Jona-Lashinio model is 
explored at flnite temperature and imaginary chemical potential with a particular emphasis on the 
confinement-deconfinement transition. We point out that the confined phase is characterized by 
a cos 3 fii/T dependence of the chiral condensate on the imaginary chemical potential while in the 
deconfined phase this dependence is given by cos Hi /T and accompanied by a cusp structure induced 
by the Z{3) transition. We demonstrate that the phase structure of the model strongly depends 
on the choice of the Polyakov loop potential U. Furthermore, we find that by changing the four 
fermion coupling constant Gs , the location of the critical endpoint of the deconfinement transition 
can be moved into the real chemical potential region. We propose a new parameter characterizing 
the confinement-deconfinement transition. 

PACS numbers: ll.30.Rd, 12.38.Aw, 12.39.Fe, 25.75.Nq 
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I. INTRODUCTION 

The exploration of the phase diagram of strongly in- 
teracting matter has received a lot of attention in recent 
years. First principle calculations of the phase struc- 
ture from the Lagrangian of Quantum Chromodynam- 
ics (QCD) is intrinsically difficult owing to the strongly 
coupled nature of the theory at large distances. Lattice 
Gause Theory (LGT) calculations provide a unique and 
powerful tool for studying QCD in the non-perturbative 
regime. Increasing computer power has recently made 
LGT simulations at almost physical quark masses possi- 
ble [l|, 0] . The results indicate that the transition from 
the confined, chirally broken phase to the deconfined, 
chirally restored phase at T '^ 160 MeV and vanishing 
baryon chemical potential, /i^ — 0, is of the crossover 
type a. 

For non-zero net baryon density, LGT calculations suf- 
fer from the so-called "sign problem" . For finite quark 
chemical potential (and Nc = 3), the statistical weight 
of the Monte-Carlo simulation becomes non positive def- 
inite due to the complex fermion determinant. This issue 
has impeded the progress in LGT calculations at finite 
densities. 

There have been several attempts to bypass the sign 
problem (J|. One interesting approach involves using an 
imaginary chemical potential, for which the fermion de- 
terminant is real and, therefore, systematic LGT simu- 
lations are possible [5|. There are two major ways for 
extracting information on the real phase diagram from 
calculations at imaginary chemical potential. One is to 
project the grand partition function Zq computed at 
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imaginary chemical potential onto the canonical parti- 
tion function 

27r 


(1) 
In spite of the difficulties involved in the evaluation of the 
oscillatory integral, there are lattice calculations aimed at 
studying the phase diagram in the temperature-number 
density plane by means of this approach [6|, |7[. An 
alternative way involves an analytic continuation from 
imaginary to real values of the chemical potential. This 
method has proven quite powerful for determining the 
critical line at fiq < ttT/S Q and this approach has been 
applied in LGT calculations |9l-[l8j. in resummed pertur- 
bation theory [1^ , as well as in quasiparticle models [20] . 
Of course, the analytic continuation requires knowledge 
of the analytic structure of the thermodynamic functions. 
Therefore effective models that share the symmetries of 
QCD are useful for testing such approaches, because a 
result obtained by analytic continuation from imaginary 
chemical potential can be confronted with the known so- 
lution at real chemical potential. 

A remarkable feature of QCD at imaginary chemi- 
cal potential is the Roberge- Weiss (RW) transition at 
fJ-q/T = 7r/3 -I- 27rfc/3, where k is an integer |5|. The 
RW transition involves a shift from one Z{3) sector to 
another in the deconfined phase. This transition is a 
remnant of the Z{3) symmetry of the pure gauge theory, 
which is explicitly broken in the presence of fermions of 
finite mass. Note that in this case the Polyakov loop is 
not an exact order parameter. Using perturbation the- 
ory [5| , Roberge and Weiss showed that this phase transi- 
tion is first-order. This was confirmed in subsequent lat- 
tice simulations 0. While it was expected that the RW 
transition is a signature of the deconfined phase 2l| , the 



transition line, which is parahel to the temperature axis 
at pLi/T = 7r/3, terminates at a temperature above the 
deconfinement transition temperature at vanishing chem- 
ical potential. Since the characteristics of the endpoint 
and its implications for the phase diagram at real ^ are 
still debated [la, [l3, I2a - l26| , it is interesting to explore 
the phase structure at imaginary chemical potential in 
an effective model. The aim of this paper is to charac- 
terize the phase structure especially of the confinement- 
deconfinement transition of QCD at imaginary chemical 
potential in the framework of an effective model which 
exhibits the relevant symmetries. 



In this work we use the Polyakov-loop-extended 
Nambu-Jona-Lasinio (PNJL) model [13, [ll. The NJL 
model [23, |30| describes many aspects of QCD related to 
chiral symmetry [31|. This model, however, lacks con- 
finement. On the other hand, thermal models with in- 
ternal gauge symmetry have been studied. These models 
reveal a RW transition [S^, |33 while chiral symmetry 
is not realized. The PNJL model is an effective model 
of QCD, which ameliorates some of the shortcomings of 
the NJL model by introducing a coupling of the quark 
field to a uniform background gauge field Aq. It has 
been demonstrated that the PNJL model reproduces the 
RW transition ^. The authors of Ref. ^ have stud- 
ied t he p hase structure of the PNJL model in detail (see 
also [2al24. l35l - l39| ). These studies indicate that various 
improvements are necessary in order to reproduce the lat- 
tice data. In this paper, applying the simplest interaction 
term in the PNJL model as introduced in Ref. [SJ], we 
focus on differences in the behaviour of the order param- 
eters dependently on the parametrizations of the effec- 
tive Polyakov loop potential. We characterize the phase 
structure qualitatively through a systematic comparison 
of the results for different Polyakov loop potentials and 
give perspectives on the nature of the phase transitions 
at imaginary chemical potential. 



We introduce a new quantity, which characterizes the 
confinement-deconfincmcnt transition based on the char- 
acteristic dependence of the chiral condensate on the 
imaginary chemical potential. The relation of this pa- 
rameter to the so called dual order parameter [40] is dis- 
cussed. 



The paper is organized as follows: in the next section 
we bricfiy review the basic properties of the QCD par- 
tition function which are relevant for this study. The 
model is introduced in Sec. lIIII and results of the numer- 
ical calculation are presented in Sec. IIVI In Sec. |Vl wc 
discuss the parameters characterizing the confinement- 
deconfinement transition and finally in Section I VII we 
summarize. 
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FIG. L Schematic phase diagram on T — 6^ plane. The 
solid line denotes the first-order Roberge- Weiss transition line 
which terminates at T = Te, while the dashed line shows the 
transition line of the chiral and confinement-deconfinement 
transitions. The arrows labeled A-C indicate paths probed 
in Sec. HVl 



II. GENERAL PROPERTIES OF THE QCD 

PARTITION FUNCTION AT IMAGINARY 

CHEMICAL POTENTIAL 

The partition function of the SU{Nc) gauge theory 
with fermions, is characterized by the number oper- 
ator N = J(Pxq^q, at imaginary chemical potential 
= M//T-/?M/. 



ZG{T,V,e) = Tr 



-IJH+iON 



(2) 



The partition function can be expressed in terms of the 
functional integral 



ZaiT^Vl 



VqVqVAf.e-'^'''^^''^'^ (3) 



where Se is the Euclidean action 



Se^ 



f dT j d?xq{-i-D-m)q--G'^ ^i^q\. (4) 
Jo J 4/5 



In ^ , the gauge and quark field obey periodic and anti- 
periodic boundary conditions in the temporal interval 
[0,/3], respectively. The particle-antiparticle symmetry 
implies, that Zq is an even function of 6, Zg{T, V, —9) = 

ZG{T,v,e). 

By performing the change of variables q -^ e^^'^^^q, 
the explicit dependence imaginary chemical potential in 
the action can be removed and converted into a modified 
boundary condition [3, [2l| 

q{x,(3)^-e'U{x,0). (5) 

Then Z(Nc) transformation 

q^Uq, A^-^UA^U-'-ii/g)id^U)U-' (6) 



with U{x,l3) = exp{2Trik/Nc)U{x,0) leaves the action 
and the functional measure invariant, but modifies the 
boundary condition 



III. PNJL MODEL 



Formulation 



q{x,P)^-e^'e^\{x,0), 



(7) 



where k is an integer [a, [21|. A comparison of Eqs. ([7]) 
and ([5]) reveals that the partition function is periodic 
with respect to finite shifts of 6 



ZG{T,V,e) = ZGiT,V, 



27rfc 



(8) 



This periodicity is a remnant of the Z{Nc) symmetry, 
the center symmetry of the SU{Nc) gauge group. In the 
presence of fermions the Z{Nc) symmetry is explicitly 
broken. In other words, the partition function is invariant 
under the transformation 6 — > 9 + 2'Kk/Nc combined with 
the Z{Nc) transformation. This symmetry was dubbed 
"extended Z{Nc) symmetry" in Ref. [3J]. It is easy to 
see that the thermodynamic potential fl = —TIuZq and 
the chiral condensate a = {qq) = —y-^^, with mg being 
the current quark mass, have the same periodicity. 

Roberge and Weiss noted the existence of the first- 
order transition a.t 9 = tt/Nc + 2nk/Nc in the decon- 
fined phase [a]. It was expected that the RW transition 
takes place at the same temperature as the confinement- 
deconfinement transition temperature Td- Lattice sim- 
ulations, however, showed that the endpoint of the RW 
transition is at a temperature Te, which is higher than 
Td- A schematic phase diagram for imaginary chemical 
potantial is shown in Fig. [T] The nature of the chiral 
and confinement-deconfinement transition at finite 9 is 
not fully understood yet. First, there is no a priori rea- 
son that these two transitions coincide. LGT simulations, 
however, show that the two transitions take place at ap- 
proximately the same temperature [1, [ifl, 13 ■ Second, 
the order of the transition depends on the quark mass and 
the number of flavors. For Nf = 2 and A'c = 3 (this case 
will be explored in this paper) Ref. [i^| shows that for 
m-^jnip ~ 0.9 the chiral and confinement-deconfinement 
transition at = 0.92(7r/3) is of the crossover type, im- 
plying that the transition at the RW endpoint is second 
order. On the other hand, Refs. [l3, [231 have demon- 
strated that for light and heavy quark masses, the phase 
transition at Te is first-order. In this case, the first-order 
RW line continues along the dashed lines and terminates 
at a critical endpoint which is located at a value of the 
imaginary chemical potential < < tt/Nc- 

In this paper, we do not attempt to obtain a fit of 
the PNJL model to lattice results at imaginary chemi- 
cal potential. Rather, we focus on exploring qualitative 
features of the phase structure of the model and their 
origin. 



Studies of the phase diagram of strongly at imagi- 
nary chemical potential within effective models, requires 
a model with the same symmetry structure as QCD. The 
most important symmetries which must be accounted for 
are the extended ^(3) and the chiral symmetries. In this 
article we employ the Polyakov loop extended Nambu- 
Jona-Lasinio (PNJL) model with Nf = 2 and A'c = 3 
[23,[2g|, which respects the above mentioned symmetries. 

The Lagrangian of the two-fiavor PNJL model with a 
four-quark interaction is given by 



C = q{i-i^,D^'-ma)q + Gs[{qqf 

-W($[yl],$*[A];r). 



{qilbfqf 



(9) 



In the covariant derivative D^^ = 9^ — iA^^ , only the tem- 
poral component Aq of the gluon field A^ = gA'^X"- /2 is 
included. The gluon field is treated as a classical back- 
ground field, whose dynamics is encoded in the effective 
potential U. The gluon field is expressed in terms of the 
traced Polyakov loop and its conjugate 



where the trace is taken over color space and 



(10) 



L{x) = Vexp 



'. / dTA^{x,T) 
Jo 



(11) 



Here A4 = iA^ and and V denotes the path or- 
dering in Euclidean time r. In the Polyakov gauge 
the matrix L reduces to the diagonal form L = 
diag(e*'^^e*'^^e-*('*l+*2)) ^. Note that in general $* 
is not the complex conjugate of $. At real chemical po- 
tential both <& and <&* are real and at // 7^ their values 
differ [4l|. At imaginary chemical potential, $* is the 
complex conjugate of $ and they have non-zero imagi- 
nary parts [33|. Therefore, in the discussion below, we 
use the notation 






me-'' 



(12) 
(13) 



for imaginary fj,. The extended Z(3) symmetry leads to 
the following properties of the Polyakov loop [3^] 



$ 9 



2TTk\ 



mo)i 



m-9)\ = m9)i 

4>{-9) = 



(14) 
(15) 
(16) 
(17) 



TABLE I. Parameters in the polynomial potential ([191 



mented l27l. |45 



To[MeV] 


ao 


ai 


a2 


as 


b3 


b4 


270 


6.75 


-1.95 


2.625 


-7.44 


0.75 


7.5 


TABLE II 


Parameters in 


the logarithmic potential (|2l[). 


To[MeV] 


ao 




ai 


^2 




bs 


270 


3.51 




-2.47 


15.22 




-1.75 



In the mean-field approximation the thermodynamic 
potential can be written in terms of quark quasiparti- 
cles with a dynamical mass M, momentum p = |p|, and 
energy Ep ^ ^Jf^TW [H [H 

n{T,V,e)^~W j -^[^{Ep-El) 



1 



ln[l + 3($ + $*e"'^^'' )e 



-^i?; 



-^PE- 



/3 
1 



^ ln[l + 3(<i>* + <^e-^''t)e-^K + g-^/^^p" 



(18) 



The term in the first line represents the ultraviolet di- 
vergent vacuum fluctuations. As usual, we introduce 
a three-momentum cutoff A in the momentum integra- 
tion to regularize the vacuum contribution. We subtract 
the vacuum term with the single-quark energy E^ — 



vw 



so that the vacuum contribution vanishes 



when the chiral symmetry is restored, as in Ref. [41|. We 
use a short-hand notation, where the imaginary chemical 
potential is subsumed in E^ = Ep ± iO / f3. The dynami- 
cal mass M is related to the current quark mass and the 
chiral condensate a = (qq) by M = mg — 2GsCr. The 
term GgC^ in the last line is due to the meson potential 
in the Lagrangian ([9]). 

So far two types of the Polyakov loop effective potential 
U have been widely used. The polynomial potental has 
as a general Z{3) symmetric form ^28. ,421 : 



W, 



poly 



b2{T) 



rp4 



$*$- 



ft L 



($*)3 



-|(**1>)' (19) 



with 



b2{T) = ao + aiy + ^2 f ^ 



, ^0 
as I Y 



(20) 



The coefficients are determined by fitting the equation 
of state and the expectation value of the Polyakov loop 
to lattice data of pure gauge theory [43, |4J| in Ref. [2^ . 
In the other widely used variant, a logarithmic poten- 
tial motivated by the strong coupling expansion is imple- 






where 



aiT) 



$*$ + 6(T) 



log{l-6$*$ + 4[$3 + ($*)3] -3($*$)2} 

(21) 



a{T) = ao + ai— + 



02 



b{T) = h 



(22) 

The above parameterization for the temperature depen- 
dency was introduced in [45] and the constants are de- 
termined by fitting lattice data of pure SU(3) theory. In 
Ref. [23] , a similar functional form in $ but different pa- 
rameterization of the temperature dependence was intro- 
duced. In this potential, however, one of the parameters 
is fixed to reproduce the simultaneous crossover transi- 
tion for chiral and deconfinement transition rather than 
the pure gauge theory except for the transition temper- 
ature To ~ 270 MeV. We refer to ^ for discussion. 
This difference makes it difficult to perform a systematic 
comparison of effect of quarks near the deconfinement 
transition. If we re-fit the parameters to reproduce the 
pure SU(3) lattice data, we expect to have similar results 
to those from the logarithmic potential (PT|) since the tar- 
get space and the transition temperature are the same. 
In this paper, we focus on the two potentials Eqs. p9)) 
and ((2T|) which equally reproduce the Polyakov loop and 
thermodynamics as well as the first order confinement- 
deconfinement phase transition. We use the parameters 
determined in Refs. [28] and [M^. For convenience, they 
are summarized in Tables |I] and [ill 

The order parameters, chiral condensate a (or dynam- 
ical mass A'l), modulus of the Polyakov loop ]$], and the 
phase of the Polyakov loop are determined numerically 
by solving the coupled equations of motion 



dn{T,V,0;M,m,, 



0, X, = M,|$],0. (23) 



The phase diagram in the T — 9 plane of the pol yno mial 
potential model pJ)) has been studied in Refs. |34l. ISSj. 
In this model, the first-order Roberge- Weiss transition 
dX 9 = 7r/3 ± 27rfc/3, the second-order chiral transition 
in the chiral limit and the crossover one at finite quark 
mass as well as the crossover confinement-deconfinement 
transition were found. However, these features, depend 
on the choice of the Polyakov loop effective potential and 
further quark interaction terms are required to reproduce 
lattice results quantitatively [2J, [S^ . 

In this paper, we restrict ourselves to the simplest 
quark-quark interaction form, as shown in Eq. ([9]) and 
focus on behavior of the order parameters in the T — 9 
plane for the polynomial and logarithmic Polyakov loop 
potentials. 



B. Some analytic insights 

Before proceeding to the full numerical computation, 
it is useful to explore the general properties of the ther- 
modynamic potential (jl8|) analytically in a few limiting 
cases. The momentum integration in Eq. (J18p can be 
carried out analytically if we first expand the logarith- 
mic terms in the integrand in powers of e~^^^ <^ 1. We 
thus find, keeping terms up to order (e~''^p)^. 



VLc^{Gs(t'^ +U)V -VLq 



2V 



3($e 



ie 



^*e- 



dpp e 



2„-/3Sp 



+ ^ {e2^»(2$* - 3$2) + e-2»''(2$ - 3$*^)} / dpp^e^^^''- 



+ {2(1 - 9$$*) cos 36* + 9e^*''$^ + 9e'^*^$*^} 



X / dpp e 



2 -3t3Ep 



(24) 



Here, ilo is the temperature independent vacuum term. 
While it is necessary to introduce a finite cutoff for 
this term, due to the non-renormalizability of the PNJL 
model, taking A — ;> cxd in the thermal part does not affect 
the qualitative features discussed below. Thus the mo- 
mentum integration of the thermal part can be carried 
explcitily resulting in modified Bessel functions Kn- 

We first consider the low temperature limit $ = $* = 
in order to explore the effect of the Polyakov loop in 
the confined phase. In this case the Polyakov loop effec- 
tive potential, U, vanishes. Furthermore, the gap equa- 
tion for the dynamical mass M, obtained from Eq. (IMl) . 
dfl/dM = 0, reduces to 

M ~ ma + —JHM'f{A/M) - mJ/(A/mo)] 
where 



f{x) = xVx^ + 1 - ln(a: + \/x^ + !)• (26) 

Note that, for 3M/T > 1, 



Ki{2,M/T) 



ttT 
6m' 



-3M/T 



(27) 



The gap equation (^5]) implies that the 9 dependence 
of the dynamical mass is completely determined by cos 3^ 
term. Consequently, M is a periodic function of 9 
with the period 27r/3, as expected. For small tempera- 
tures, quark degrees of freedom are strongly suppressed. 
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since, for a vanishing Polyakov loop, only 



FIG. 2. Temperature dependence of the chiral condensate 
normalized by the value at T = and the Polyakov loop at 
vanishing chemical potential. Upper panel shows the result in 
the chiral limit while the lower stands for the physical quark 



three-quark clusters survive. The chiral phase transition 
takes place when the thermal contribution in the gap 
equation is of the same order as the vacuum term. If the 
condition $ = $* = is strictly enforced, the chiral tran- 
sition is shifted to very high temperatures. Because the 
thermal excitation of quarks in this limit is possible only 
in three quark clusters, the resulting thermodynamics is 
qualitatively similar to that of the nucleonic NJL model, 
which also yields a very high chiral transition tempera- 
ture fiTJ . 

In the (naive) high-temperature limit, $ = $* -^^ 1, 
the gap equation reduces to that of the ordinary NJL 
model, 

M = mo + -^[M'fiA/M) ~ mj]/(A/mo)] 

TT 

2AGsM^T COS 9 , , , 
—^ KiiMT). 28 

Now, the 9 dependence of the dynamical mass is deter- 
mined by cos 9 and the thermal factor is proportional to 
Ki (M/T) , appropriate for the thermal excitation of sin- 
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FIG. 3. Phase of Polyakov loop as a function of 6 for various 
temperatures in the chiral hmit. As in Fig. [J] each hne stands 
for the results of the logarithmic Polyakov loop potential and 
the polynomial one at the different temperatures, respectively. 



gle quarks. In this case, the dynamical mass has lost the 
original periodicity of the partition function (|8|) , which is 
respected in the low-temperature limit (j25|) . 



At 9 — -K 12 the thermal contribution in ([^51) vanishes 
and consequently the dynamical mass equals its vacu- 
um value, irrespective of temperature. Hence, in this 
approximation the phase boundary is shifted to higher 
temperatures as the imaginary chemical potential is in- 
creased, and eventually approaches T = cxd in the limit 
Q — > 7r/2. In general, the thermal contribution is of the 
form - Yjn ^ri cosneKi{nM/T) with L„ > 0. Due to the 
higher order terms, the transition temperature remains 
finite, but the positive curvature of the phase boundary 
persists |^. On the other hand, for real chemical poten- 
tials the leading thermal contribution is proportional to 
T cosh /3/i. Since cosh x is an increasing function of x, this 
implies that the chiral transition temperature decreases 
as the real chemical potential grows. 

In the high-temperature limit, as implemented above, 
the original periodicity of the partition function is lost 
because the phase of the Polyakov loop is neglected. In 
fact, at imaginary chemical potential, the high tempera- 
ture limit of the PNJL model is in general not the NJL 
model. The Roberge- Weiss transition, characterized by 
discontinuous jumps of the phase 0, preserves the peri- 
odicity 27r/3 in the deconfinement phase. 



IV. BEHAVIOR OF THE ORDER 

PARAMETERS AT IMAGINARY CHEMICAL 

POTENTIAL 

We now discuss the characteristics of the order param- 
eters obtained by solving the full gap equation Eq. ([23]) . 
Besides the Polyakov loop potentials, given in the pre- 
vious section, the model has three parameters in the 
fermion sector. We use 



Gs 
A 



5.498GeV"^ 
0.6315GeV, 



(29) 
(30) 



^ Owing to the reflection symmetry and periodicity, it is sufficient 
to consider the interval < 9 < n/3. 



which reproduce the vacuum pion mass and pion decay 
constant at zero temperature and density, when the cur- 
rent quark mass is fixed to the value ttiq = 5.5 MeV j3l| . 

In what follows, we compare the results obtained for 
a finite pion mass with those corresponding to the chiral 
limit, mo = 0. In the latter case, the model belongs to 
the universality class of the three dimensional 0(4) spin 
model and exhibits a second-order phase transition at 
finite temperature and small values of the real chemical 
potential. 

With the parameter set given above, we find the chiral 
condensate and the Polyakov loop shown in Fig. [5] The 
dependence of the order parameters on the temperature 
shows that in the chiral limit, the chiral transition is in- 
deed second-order, while the confinement-deconfinement 
transition is of the cross over type. For finite quark mass, 
mo = 5.5 MeV, shown in the lower panel, the chiral order 
parameter and the Polyakov loop both exhibit smooth 
crossover transitions. Thus, the explicit symmetry break- 
ing induces a qualitative change of the chiral condensate, 
while for the Polyakov loop this dependence is negligible. 
Furthermore, a comparison of the two parametrizations 
of the Polyakov loop effective potential shows that the 
transition is smoother for the polynomial potential than 
for the logarithmic one. 



A. Behavior across = vr/S 

We now consider the order parameter as a function of 
8 close to 9 = 7r/3. In Figure [3] we show the phase of the 
Polyakov loop (/> as a function of 6 in the chiral limit. We 
do not show results for non-zero quark mass, since the 
results are indistinguishable from those shown in Fig. |3| 

At low temperatures (arrow A in Fig. [IJ , the phase of 
the Polyakov loop changes smoothly from at = to 
— 7r/3 at 6 = 7r/3. Subsequently, the phase continues to 
decrease and finally approaches to — 27r/3 a.t 9 = 2tt/3, 
as required by the symmetry, Eq. (1161) . This behavior is 
independent of the choice of U. At temperatures beyond 
Te, the transition is first order (arrow C in Fig. [T|). For 
instance at T = 280 MeV, the phase jumps from = to 
— 27r/3 at 9 — tt/3. This is the Roberge- Weiss transition 
[5|, where the phase of the Polyakov loop jumps from one 
Z{3) sector to another. The RW transition is common 
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to both parametrizations of U. This is natural, since the 
RW transition is a consequence of of the Z{3) symmetry 
of the pure gauge theory, which is incorporated in both 
potentials. The detailed behavior around Te is, however, 
different between the two potentials. Thus, at T = 250 
MeV for the logarithmic potential, which is below the 
endpoint of the RW transition (Te ~ 255 MeV), the 
phase is discontinuous ai 9 ^ 7r/3. This implies that the 
phase boundary, which is crossed by arrow B in Fig. [T] is 
first order at this temperature. By contrast, in the poly- 
nomial case the phase is a smooth function of 9 at any 
temperature below the RW endpoint {Te — 275 MeV). 
As we discuss below, this reflects the different order of 
the RW endpoint for the two potentials. 

We note that the logarithmic potential is defined in a 
limited domain, characterized by positivity of the argu- 
ment of the logarithm, while for the polynomial potential 
there is no such restriction. In fact, for high temperatures 
(e.g. T = 280 MeV) the Polyakov loop, plotted for the 
polynomial potential as a function of 9 in the complex 
<& plane, leaves the so called target space, defined by re- 
quiring that the logarithmic potential is real j48|. 

In Figure 2] we show the modulus of the Polyakov 
loop l^l as a function of 9. At high temperatures, the 
RW transition is manifested by a cusp at 9 = 7r/3, 
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while at lower temperatures \^\ varies smoothly for 
both parametrizations of the Polyakov-loop potential. 
However, at intermediate temperatures, the two poten- 
tials yield qualitatively different results, as illustrated 
by the discontinuities in |$| obtained for the logarith- 
mic potential near 9 — 7r/3 at T = 250 MeV. For the 
polynomial potential we find a continuous confinement- 
deconfinement transition at imaginary chemical poten- 
tial, while for the logarithmic potential the transition is 
first order at intermediate temperatures. We return to 
this point in the following subsection. Here we note only 
that the first order transition is reflected also in a sudden 
change of the phase at T = 250 MeV, as shown in Fig. [31 

Within the PNJL model, the transition from one Z{3) 
sector to another can be understood in the following 
way. In the high temperature limit, the dominant con- 
tribution to the thermodynamic potential is given by the 
single quark excitation term in Eq. (I24p (the first term 
in the square bracket), which yields a contribution to 
il r^ — cos{9 + (j)). The only additional i/i-dependent term 
in f2 is the Polyakov loop potential U. At high tempera- 
ture, U has the three local minima at ^ = and ±27r/3. 
For each value of 9, the physical vacuum is obtained by 
finding the absolute minimum of the two terms. As illus- 
trated in Fig. [SI the physical vacuum changes from one 
minimum to the next as 9 crosses 7r/3 -I- 27rfc/3. While 
both potentials have the periodicity A9 — 27r, when 
is artificially fixed in one Z{3) sector, the complete ther- 
modynamic potential acquires the periodicity A9 — 27r/3 
owing to the (f> dependence of U. 

Similarly, the chiral condensate a is also continuous 
near 9 = 7r/3 at temperatures lower than Te, irrespective 
of the quark mass and the potential, as shown in Fig. [5] 
At temperatures higher than Tg;, it develops a cusp for 
both parametrizations, a manifestation of the RW transi- 
tion. Furthermore, at temperatures above T^(9 = 7r/3), 
a vanishes for any value of 9 in the chiral limit, as shown 
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for T — 340 MeV. For a finite quark mass, the cusp per- 
sists to temperatures much higher than T^ but finally 
disappears, as shown in Fig. [6] for toq = 5.5 MeV and 
T = 500 MeV. 

In the chiral limit, for temperatures in the interval be- 
tween T^{0 = 0) and Tj^{9 = 7r/3), there is a second 
order chiral transition at non-zero 9, as shown in Fig. [51 
(see also arrow B in Fig. [I}. For finite quark masses, 
this transition is of the cross-over type, as illustrated in 
the lower panel of Fig. [S] for mo = 5.5 MeV. Also the 
temperature dependence of a is clearly different for the 
two potentials. In particular, there is a discontinuity in a 
near 6 = tt/B at T = 250 MeV for the logarithmic poten- 
tial. The values of temperature and imaginary chemical 
potential corresponding to the discontinuity are identical 
to the ones obtained for the Polyakov loop. 



B. First-order phase transition at T < Te 

In this section we focus on the first order transition 
found in a limited temperature range below the Roberge- 
Weiss transition for the logarithmic potential. In Fig. [7] 
we illustrate this result at T = 250 MeV. In each panel 
two lines are shown: one is the solution of the gap equa- 
tions approaching the transition from small 0, while the 
other is obtained by approaching from the opposite side. 
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The existence of two solutions in a certain range of 9 
shows that there are two local minima in the thermody- 
namic potential. The first-order phase transition takes 
place at the value of 9, where the thermodynamic poten- 
tial in the two local minima is degenerate. At T = 250 
MeV, this happens at 6* = 0.911(7r/3). The lines ter- 
minate where the corresponding minimum disappears. 
Thus, the system exhibits hysteresis, a characteristic of 
a first-order phase transition. 

Although this transition is related to the confinement- 
deconfinement transition, the discontinuity is reflected 
also in the chiral condensate a, owing to the coupling 
between the Polyakov loop and the chiral order parame- 
ter. 
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C. The RW endpoint 

The existence of the first-order transition, discussed 
in the previous section, is closely related to the char- 
acteristics of the RW endpoint. In Fig. |S] we show the 
temperature dependence of the Polyakov loop along the 
line 9 = 7r/3 (cf. Fig. [IJ for the two potentials in the 
chiral limit. The nature of the RW endpoint, as char- 
acterized e.g. by the phase of the Polyakov loop, dif- 
fers between the two potentials. While the polynomial 
potential yields a continuous transition, the logarithmic 
one exhibits a discontinuity in the phase and magnitude 
of the Polyakov loop. 

Above the RW endpoint (T ~ Te), the phase of the 
Polyakov loop can take two values on the 9 — tt/3 line, 
corresponding to different Z(3) sectors (cf. Fig. [5]). Thus, 
in the case of the polynomial potential, the phase bifur- 
cates smoothly at the RW endpoint, while for the loga- 
rithmic potential, the phase changes discontinuously at 
this point. In the former case, the RW endpoint is a sec- 
ond order point, while in the latter it is a triple point. 
We note, however, that a different parameterization for 
the logarithmic potential [23] also yields a second order 
RW endpoint ^. 

Consequently, the characterization of the RW endpoint 
depends on the parametrization of the Polyakov loop po- 
tential, but, within the framework considered here, it is 



independent of the value for the quark mass for the log- 
arithmic potential. On the other hand, in LGT calcu- 
lations it is found that the nature of the RW endpoint 
depends crucially on the quark mass; for both two and 
three flavor QCD it is a second order endpoint at inter- 
mediate quark masses and a trip le point for large and 
small masses [ij, 123 ■ In Ref . [2J] it is argued that in or- 
der to reproduce the quark mass dependence of the RW 
transition, a $ dependent fermion coupling, motivated 
by functional renormalization group analyses [4^, [50] , is 
required. This indicates that in QCD the interplay be- 
tween chiral symmetry breaking and confinement is more 
complicated than in the present model. 

Finally, in Fig.|9]we summarize the results on the phase 
diagram in the T — 9 plane for the two Polyakov loop 
potentials. The (pseudo-)critical temperature for the de- 
confinement transition corresponds to a maximum of the 
temperature derivative of the modulus of the Polyakov 
loop d\^\/dT. The critical end point (CEP), obtained 
for the logarithmic potential at < 6* < 7r/3, is a conse- 
quence of the triple point at T = Te and 9 = tt/3. 

In comparison of the two potentials, we have seen that 
the phase transitions at imaginary chemical potential are 
shifted to higher temperatures compared to those at real 
chemical potential. This implies that dynamical quark 
mass becomes heavier at fixed temperature (see Fig. [5]) 
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thus the Polyakov loop potential W, which is indepen- 
dent of 9, makes a dominant contribution to the ther- 
modynamic potential. Furthermore, at imaginary chem- 
ical potential, the target space of the Polyakov loop is 
probed through the change of the phase (/). Therefore, a 
comparison of the resulting phase diagram at the imagi- 
nary chemical potential region with that obtained in LGT 
calculations, yields important constraints on the effective 
Polyakov loop potential. 



D. Critical endpoint of confinement-deconfinement 
transition 

In Ref. fSl'l it was found that as the pion mass in the 
quark-meson model is reduced from its physical value, an 
additional critical endpoint appears on the phase bound- 
ary at small (real) values of the chemical potential. Since 
the coupling to the Polyakov loop is not accounted for 
in [511, the additional CEP is associated with the chiral 
phase transition [j. In this section we explore the de- 
pendence of the confinement-deconfinement CEP, which 
appears at imaginary /i for the logarithmic potential, on 
the model parameters. We find that the location of this 
CEP depends on the four fermion coupling constant Gg- 

In Fig. [lU] we show the phase diagram of the model 
in the chiral limit for different values of Gg- We in- 
clude both real and imaginary values of /i by showing 
the phase boundaries in T — /i^ plane. The upper pan- 
els show cases where Gs is smaller, while the lower ones 
show cases where it is larger than the reference value 
(P5)l . Lines appearing for /x^ < — (ttT/S)^ (the boundary 
is indicated by the dotted line) are images of those in the 
region — (ttT/S)^ < M^ < 0; the mapping is defined by 
the periodicity of the partition function. 

We identify the (pseudo)critical temperature by find- 
ing the the maximum of the derivative of the correspond- 
ing order parameter with respect to temperature. For 
real values of the chemical potential, the Polyakov loop 
and its conjugate are real but take on different values 
4l|. Here we use d^/dT for the definition of the de- 
confinement transition. A different definition, based e.g. 
on the Polyakov loop susceptibility [4l| , would lead to a 
slightly different value of the pseudocritical temperature 
of the crossover transition. 

For the first order transition at large /i^ one finds dou- 
ble peaks in d^/dT (cf. Fig. [TU]) . Here we identify 
the position of the transition with the maximum which 
smoothly extrapolates to the deconfinement transition at 
vanishing chemical potential and to the chiral transition 
at T = 0, where the peak structure is simple. We note 



^ The fact that the model calculation of Ref. [Slll yields a first order 
chiral transition at both small and large values of fi in the chiral 
limit, is presumably due to the neglect of the fermion vacuum 
loop [12. 



that any ambiguity in the location of the phase boundary 
does not affect the discussion below. 

Qualitatively the features of the phase diagram can 
be classified as follows. The NJL sector has a critical 
coupling G" = ■k'^ l{2NrA^) for the gap equation to have 
a nontrivial solution |3l| . This implies that, with the 
present three-momentum cutoff, there is no spontaneous 
breakdown of chiral symmetry for Gs < 4.125 GeV^^. 
Therefore, in the upper-left panel {Gs — 3.5 GeV^^), 
the system is everywhere in the chirally symmetric phase. 
In this case the RW endpoint is still a triple point, and 
the CEP of the deconfinement transition is close by, at 
9 = 0.957r/3. 

For Gs > G", there is a chiral transition at vanishing 
chemical potential. As seen in the upper-center panel, 
there is a precursor at imaginary chemical potential for 
Gs slightly smaller than G"'. The chiral symmetry is 
spontaneously broken in a small region at intermediate 
temperatures adjoining the /Lt^ = —{ttT/S)'^ line. This 
behavior can be understood along the lines presented 
in Sec. IIIIBI Although the Boltzmann approximation, 
Eq. (I24p . might not be a good approximation since the 
system is in the chirally symmetric phase even at low 
temperatures, the Polyakov loop is small so that the ther- 
mal contribution is dominated by the cos 39 term as in 
Eq. (|25|) . Since this term is positive in the 7r/6 < 9 < 7r/2, 
it adds to the vacuum term and a non-trivial solution ap- 
pears at finite temperature where the system enters the 
broken phase. 

As the temperature is increased further, however, the 
Polyakov loop is non-zero and the one- and two- quark ex- 
citations contribute to the gap equation, driving the sys- 
tem back into the symmetric phase. Consequently, near 
the RW endpoint, the chiral and deconfinement transi- 
tions occur simultaneously and the chiral transition is 
also of first order. Note that the lower endpoint of the 
chiral transition follows the ^^ = —{irT/i)'^ line and ar- 
rives at the origin when Gs = Gf. For Gs beyond this 
value, the chiral transition line enters the /z^ > half- 
plane. 

As Gs is increased beyond Gs = Gg'^, the location of 
the deconfinement CEP moves to larger /z^ and finally 
reaches real values of the chemical potential at G^ — 6.5 
GeV~^. At the same time, the chiral transition line 
moves to larger T and /i^. The behavior of the CEP 
can be directly related to changes of the chiral transi- 
tion with increasing Gs- Increasing Gs leads to a larger 
dynamical mass M = mo — 2GsCr in the chirally broken 
phase. Moreover, \a\ increases with Gs since a stronger 
scalar coupling leads to larger quark condensate. This 
raises both the dynamical mass and the chiral transition 
temperature. On the other hand, a modified Gs does 
not alter the Polyakov loop sector of the model. Conse- 
quently, near the deconfinement transition, the dynami- 
cal mass of the quarks increase with Gs and the system 
approaches a pure gauge theory owing to the thermal 
suppression of quark degrees of freedom. Thus, as shown 
in Fig. 1111 the first order phase transition of the pure 
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gauge theory is recovered at Gg = 6.3 GeV ^, only 15% 



above the reference value (|29|) . 

One may wonder why this mechanism is not effective 
for the polynomial potential, since it also exhibits a first 
order phase transition in the absence of fermions. In 
fact, we find that the polynomial potential does show the 
same behavior, but at much larger values of the scalar 
coupling. The RW endpoint, which is second order at 
Gs = 5.498 GeV~^, is first order transition starting at 
Gs = 12.4 GeV^^ and the deconfinement CEP reaches 
n'^ = at Gs = 25 GeV~^, 4.5 times larger than the 
reference value. 



The origin of this quantitative difference between the 
polynomial and logarithmic potentials is due to the much 
weaker first order phase transition (smaller discontinuity 
in $) exhibited by the polynomial potential in the heavy 
quark limit. At Gs = 25 GeV~^ we find a dynamical 
quark mass of about 2.5 GeV. Thus, for a quark mass less 
than 2.5 GeV the first order transition is smoothened to a 
crossover transition. By contrast, for the logarithmic po- 
tential this happens at a much smaller dynamical quark 
mass of 0.4 GeV, owing to the much stronger underlying 
first order transition. 

We note that a first order confinement-deconfinement 
transition emerges at real chemical potential also in the 
large Nc limit of the PN JL model |53l , as explored in the 
context of the recently proposed quarkyonic phase j54| . 
Indeed, the effect of strengthening Gg is similar to that of 
increasing Nc since Gs and Nc appear in the factor GsNc 
in the gap equation for the dynamical mass (see Eq.(16) 
in Ref. l53|). While we suppress the quark contribution 
to the thermodynamics by increasing the dynamical mass 
by means of G^, a large value of Nc is accompanied by 
a 1/Nc suppression of the quark contribution. Both pro- 
cedures yield a giuon dominated system and thus give a 
first order confinement-deconfinement transition. 

Note, however, that the two procedures differ in de- 
tail. Increasing Gs at fixed Nc preserves the Polyakov 
loop potential but modifies the quark mass, while in the 
large Nc limit at fixed GgNc the quark mass remains 
unchanged but the Polyakov loop potential is modified. 
This means that when we increase Gs at fixed Nc we 
change the characteristic scale of the chiral symmetry 
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FIG. 12. Chiral condensate a = (gg) as a function of the 
twisted angle ip. 



breaking. Although this does not correspond to the phys- 
ical situation, since QCD has unique scale, Aqcd, our 
result could be useful for exploring the interplay between 
the chiral phase transition and deconfinement. 



V. DUAL PARAMETER FOR THE 

CONFINEMENT-DECONFINEMENT 

TRANSITION 

In this section, we consider dual parameters which 
capture the characteristic feature of the confinement- 
deconfinement transition discussed above. Recently a 
dual parameter has been introduced by considering a gen- 
eralized boundary condition for fermions 



q(x,^) = e^q(x,Q). 



(31) 



Here iy9 is so-called the twisted angle. The dual quark 
condensate S^") is defined as the n-th Fourier coefficient 
of the chiral condensate as a function of the twisted angle 



S(")(T) = - / " ^e-™^a(r,^). 
Jo 27r 



(32) 



The chiral condensate (T(r, 1^9) is defined in terms of kp by 

m 



cr(^) 



(Tr[(mo + i5^)-i]) 



G- 



(33) 



In Ref. [401 ■ this quantity was introduced based on the 
lattice regularization. The ^p dependence of Tr[(TOo + 
D^)~^\ can be written down explicitly by using the link 
variable. It reduces to the ordinary chiral order parame- 
ter in the limit of m ^ and Y ^ 00. An implementa- 
tion in the continuum theory has been done in the frame- 
work of Dyson-Schwinger equation |55| . The most inter- 
esting quantity is that of n = 1, which is called dressed 
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Polyakov loop. Because of a relation to the ordinary 
Polyakov loop, it can be regarded as an order parameter 
of the confinement-deconfinement transition. 

From the fermionic boundary conditions pip , one im- 
mediately finds that this is equivalent to introducing the 
imaginary chemical potential [Eq. ([5])]. The only dif- 
ference is that (/9 — ±7r corresponds to the usual anti- 
periodic boundary condition in the twisted angle while 
= does so in the imaginary chemical potential. In 
this case, the relation between (p and d is given by just a 
shift of ±7r. 



(34) 



Lp = B±-K (mod 27r/iVc). 



Furthermore, one notes that the dual condensate ([32]) is 
quite similar to the canonical partition function ([T]) . 

However, the LGT calculations demonstrate that oiip) 
exhibits quite different behavior from that for the imag- 
inary chemical potential |56j . cr((^) shows a periodicity 
of 27r in (/?, not 27r/3 which is required by the RW peri- 
odicity. The origin of this difference is the expectation 
value of the operator (• • • )g in Eq. ([55]). The subscript 
G denotes the path integral over the gauge field with the 
fermion determinant which follows the ordinary bound- 
ary condition. The change of the boundary condition 
([5T|) applies only to the Dirac operator. In the case of the 
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FIG. 15. Same as Fig. [131 but for mo = 5.5 MeV. 



imaginary chemical potential, different 9 gives a different 
fermion determinant while in the case of the twisted angle 
the background field does not change with the boundary 
condition. Therefore ip dependence of the chiral con- 
densate differs from 9 dependence. In a PNJL model, 
the authors of Ref. [571] use the value of the Polyakov 
loop calculated at = to obtain the chiral condensate 
(t((^). This prescription corresponds to varying only the 
fermionic boundary condition without changing gluonic 
background. We will follow the same prescription below. 
Since the periodicity 27r/3 in the imaginary chemical po- 
tential is preserved by the RW transition, which is an ef- 
fect of the Polyakov loop in the context of PNJL model, 
the relation p4l) holds for the normal NJL model calcu- 
lation which does not couple to ^(3) field. In spite of the 
absence of confinement in the NJL model, one sees that 
behavior of the dual chiral condensate is quite similar to 
one obtained from lattice QCD and PNJL model [58| . 

Figure [TH shows the chiral condensate as a function 
of the twisted angle (/s, obtained by the same method as 
used in Ref. [53|. The periodicity is no longer 27r/3. One 
also notes that there is broken phase at region far from 
Lp = -K even at high temperatures. This is in contrast 
to the case of the imaginary chemical potential shown 
in Fig. |6]and similar to what was expected from the gap 
equation of the NJL model, Eq. (pS)) . Indeed, at </? = 7r/2 
and 37r/2, which correspond to 9 = T7r/2, the thermal 
term vanishes in Eq. (I28|) . resulting in the almost temper- 



ature indepedent chiral condensate. On the other hand, 
there are no qualitative difference between the logarith- 
mic potential and the polynomial one. The reason is that 
the Polyakov loop enters in the gap equation only as a 
constant determined at <p = tt at each temperatures. 

Let us introduce a new dual parameter by using 9 in- 
stead of tp such that it captures the characteristics in the 
9 space. We define 



ZTT 



t/3 



d9e 



V(r,. 



(35) 



-ir/3 



The integration region is changed to [— 7r/3 : 7r/3] with re- 
spect to the periodicity 27r/3. As discussed in Sec. lIIIBl 
the physical meaning of the periodicity is different in con- 
finement and deconfincmcnt phase. In the confinement 
phase, periodicity 27r/3 is coming from cos 30 which char- 
acterizes the confinement of the quarks. On the other 
hand, deconfinemcnt phase is characterized by cos 9 with 
discontinuity at = 7r/3 + 27rfc/3 caused by Z{i) tran- 
sition which preserves the periodicity 27r/3. Therefore, 
we expect, that S^ and E^ demonstrate characteristic 
behavior of the confinement-deconfinement transition. 

In Figs. [13] and [M] we compare the three kinds 
of the characteristic parameters of the confinement- 
deconfinement transition and their derivatives with re- 
spect to temperature as functions of temperature. We 
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FIG. 16. Same as Fig. [H but for mo = 5.5 MeV. 



normalized the dressed Polyakov loop and the modified 
dual parameter ()35p so that they tend to at low temper- 
ature and to unity at high temperature by I]norm(T) = 
[E(T) - I](Ti)]/[E(T2) - S](Ti)] where we used Ti = 50 
MeV and T2 = 1 GeVO The same normalization con- 

in) 

slants are applied to the derivatives. Note that S^ be- 
fore normalization vanishes at temperature higher than 
the chiral transition temperature at = 7r/3 in the chiral 
limit since (7(T, G) does so, as seen in Fig. [Bl After the 
normalization, it smoothly approaches to unity as one 
sees in Figs. USHIl 

One observes qualitatively that both dual parameters 
show similar behavior to the Polyakov loop. The dressed 
Polyakov loop is almost parallel to the Polyakov loop in 
the temperature region covered in the figure. On the 
other hand, the modified dual parameter Eg reaches 
quickly its limiting value. The behavior around the tran- 
sition temperatures reflect the difi^erence between two po- 
tentials. In the case of the logarithmic potential, one 
sees different structures for each derivative of the order 
parameters. While the derivative of the Polyakov loop 
d\^\l dT shows the only peak corresponding to the pseud- 
ocritical temperature, the dressed Polyakov loop dY} /dT 
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FIG. 17. Dual order parameters of n = 3. Upper panel shows 
the result from Eq. p2|) while lower one displays that from 
Eq. dSH). 



Note that Eg (T = 0) does not vanish since the integration is 
from — 7r/3 to it/3, not from —tz to n. 



exhibits the two peak structure. One agrees with that 
of the d^/dT and the other corresponds to the chiral 
transition (see Fig. [2]). The derivative of the modified 
dual parameter dT,^ /dT also shows a peak for the chi- 
ral transition. However, the deconfinement appears only 
as a shoulder. The polynomial potential shows a broader 
peak in d^/dT reflecting the weaker nature of transition. 
However, maxima in dual order parameters associated 
with the crossover transition do not appear. 

Inclusion of small quark mass, TOq = 5.5 MeV slightly 
modifies the behavior of the characteristic parameters, as 
expected from the difference in a. FiguresfTSlandfTBlshow 
the three parameters and their derivatives for tuq = 5.5 
MeV. While at finite quark mass there is little differ- 
ence in the behavior of the parameters seen in toq = 
case, distinct peak structures show up in their deriva- 
tives. For the logarithmic potential, the peak associated 
with the chiral transition seen in the chiral limit does not 
exist in both dual parameters. The remnant of the chiral 
transition appears only in the modified dual parameter 
as a dip. On the other hand, polynomial potential case 
exhibits much broader peak, which corresponds to de- 
confinement in d\^\/dT and to chiral transition in dual 
parameters. 

Figure [17] shows the dual parameter ([32]) and (|35l) at 
n = 3. In all the cases, one sees that they start to in- 
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crease, exhibit peak structure and then decreasing. This 
is common for both E^^^ and Eg . Such behavior of Eg 
can be understood by analysing cr(T^ 9) which is shown 
in Fig. ini At low T, it oscillates according to cos 39 with 

(3) 

an amplitude given by the thermal factor. Then Eg , 
which is proportional to the amplitude, increases with T. 
AtT >T^, a{9, T) vanishes at 6* = [0 : 9^] thus only the 
region [9c ■ 7r/3] contributes to the integral. Therefore, 
Eg has a peak aX T — T^ and then decreases. This 
does not relate to the deconfinement phenomenon and is 
common for E^"^) and Eg ' . 

If one focuses on the difference between the logarith- 
mic an the polynomial potential, however, one finds 
a remnant of the deconfinement-confinement transition 
in the behavior of Eg . As we have seen in Fig. [71 
cr(T, 9) has a discontinuity induced by the first-order 
confinement-deconfinement transition. It is reflected to 
the non-monotonic behavior of Eg ' between T = 0.24 
and T = 0.26 GeV in the case of logarithmic potential. 
Indeed, the shoulder seen at T = 0.24 GeV corresponds 
to the temperature of CEP, at which (j{9,T) has a discon- 
tinuity, and the second inflection point reflects the RW 
endpoint which is the triple point in this case. This indi- 

(3) 

cates that the newly introduced dual parameter Eg has 
sensitivity to the confinement-deconflnement transition 
at imaginary chemical potential. 



VI. SUMMARY 

We have studied the confinement-deconfinement tran- 
sition in the PNJL model at imaginary chemical potential 
with the simplest interaction. We discussed the origin 
of the characteristic periodicity 27r/3 of the order pa- 
rameters. It is characterized by cos 30 in the confined 
phase while it is due to cos 9 with the RW transition at 
9 = 7r/3 -I- 27rfc/3 induced by the change of the phase 
of the Polyakov loop. We also explored the results from 
different Polyakov loop potentials. We found that the 
property of the confinement-deconfinement transition de- 
pends on the choice of the potential in spite of the fact 
that both potentials exhibit the first-order phase transi- 
tion in the absence of quarks. Remarkable differences are 



seen in both the RW endpoint and the behavior of the 
Polyakov loop at finite 9. In the case of the logarithmic 
potential which has more relevant domain, we find that 
the confinement-deconfinement transition becomes first 
order near 9 ~ 7r/3 and there is a critical endpoint of the 
transition at imaginary chemical potential. We also find 
that the location of the CEP moves with the four fermion 
coupling Gs and it reaches real chemical potential region 
by increasing Gg- This behavior can be understood by 
the suppression of the quark contribution since increas- 
ing Gs implies larger (qq) which quantifies the dynamical 
quark mass. At large coupling, the existence of the CEP 
is independent of the choice of the potential. However, 
the polynomial potential requires larger Gs because it ex- 
hibits much weaker first order transition. Consequently, 
it seems that the order of the deconfinement transition 
is determined by the size of the gap A<I> in the Polyakov 
loop potential and the quark condensate (qq) . 

The first order phase transition influences the behavior 
of the chiral condensate as a sudden jump at the criti- 
cal imaginary chemical potential. We proposed modi- 
fied dual parameters using the imaginary chemical po- 
tential based on the analogy to the twisted angle in the 
dual order parameters. Comparing the n = 1 case with 
the Polyakov loop and the dressed Polyakov loop, we 
found that each parameters has different sensitivity to 
the phase transitions. We showed that n — 3 has a char- 
acteristic behavior owing to the first order confinement- 
deconfinement transition at intermediate 9. We expect 
relevance of our study in understanding the QCD phase 
diagram. 
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